The IRIS Far-Infrared Galaxy Survey : 
Expected Number Count, Redshift, and Perspective 



TSUTOMU T. TAKEUCHI 12 , HIROYUKI HIRASHITA 12 , KOUJI OHTA 1 , 
TAKASHI G. HATTORI 1 , TAKAKO T. ISHII 13 , 

1 Department of Astronomy, Faculty of Science, Kyoto University, Sakyo-ku, Kyoto 
^ i 606-8502, JAPAN 

2 Research Fellow of the Japan Society for the Promotion of Science 
3 Kwasan and Hida Observatories, Yamashina-ku, Kyoto University, Kyoto, 

Q 607-8471, JAPAN 
Electronic mail: takeuchi, hirasita, ohta, hattori, ishii@kusastro.kyoto-u.ac.jp 

AND 

HIROSHI SHIBAI 4 

CN . Division of Particle and Astrophysical Sciences, School of Science, Nagoya University, 

(N 
(N 



Chikusa-ku, Nagoya, 464-8602, JAPAN 
Electronic mail: shibai@toyo.phys.nagoya-u.ac.jp 



OO 
ON 

' ABSTRACT 

6 

+l5 . Infrared Imaging Surveyor (IRIS) is a satellite which will be launched in 

the beginning of 2003. One of the main purposes of the IRIS mission is an 
all-sky survey at far-infrared (FIR) with a flux limit much deeper than that 
of IRAS. In order to examine the performance of the survey, we estimated 
the FIR galaxy counts in four (50, 70, 120, and 150 /im) bands based on 
some models. We adopted a multicomponent model which consists of cirrus 
and starburst components for galaxy spectra, and the nearby FIR luminosity 
function derived from that of IRAS galaxies. We derived the number counts, 
redshift distributions, and infrared diffuse background radiation spectra for i) no 
evolution, ii) pure luminosity evolution, iii) pure density evolution with q = 0.1 
and 0.5. We found that a large number of galaxies ( ~ a few x 10 6 in the whole 
sky) will be detected in this survey. With the aid of a vast number of detection, 
we will detect the effect of galaxy evolution, and evaluate the amplitude of 
evolution at least in the nearby universe in the IRIS survey, though it will be 
still difficult to constrain which type of evolution takes place from the number 
count alone. We also studied the estimation of redshifts of detected galaxies 
by their infrared colors alone. Although significant contamination takes place 
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among nearby faint galaxies and high-z ones, we found that rough estimation of 
galaxy redshift can be practicable by jointly using present and future optical 
surveys. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: 
starburst — infrared radiation 



1. INTRODUCTION 

One of the important problems in astrophysics is to trace back the galaxy evolution 
to the epoch at which galaxies formed. Recently Steidel and collaborators have pioneered 
and widely utilized a multicolor technique to select young galaxies by use of Lyman-break 
dropout, and discovered a class of actively star-forming galaxies in high-redshift (z J> 3) 
Universe (Steidel, Pettini,& Hamilton 1995; Steidel et al. 1996; Lowenthal et al. 1997). 
Some lens-magnified normal galaxies at high redshift also show rather high star formation 
rate (Ebbels et al. 1996; Yee et al. 1996; Franx et al. 1997; Trager et al. 1997). 
Furthermore, Ohta et al. (1996) and Omont et al. (1996) detected CO emission lines in 
the quasar BR1202-0725 at redshift z = 4.69, suggesting the existence of a large amount of 
molecular gas and dust in the object. All of these results imply that strong star formation 
and synthesis of heavy elements have already occurred in such an early stage of the 
Universe. The violent initial bursts of star formation must have produced large amount 
of dust, so that the optical deep surveys of protogalaxies inevitably suffer from the bias 
induced by severe dust extinction. The dust particles reprocess the starlight, and bulk of 
the bolometric luminosity of these young galaxies will be emitted at far-infrared (FIR) - 
submillimeter (sub-mm) wavelengths. 

The FIR emission from galaxies, especially from those are so faint that we are unable 
to resolve, produces the cosmic infrared background radiation (CIRB). Therefore the CIRB 
provides important information on the past star formation history of galaxies including the 
inaccessible sources (e.g. Bond, Carr, & Hogan 1986; Burigana et al. 1997). The recently 
reported detection of an isotropic diffuse light at A ~ 100 — 200/im has given this field new 
impetus (Puget et al. 1996; see also Fixsen et al. 1998; Hauser et al. 1998). 

Meanwhile, galaxy evolution at low reshift (z <> 1) has also crucial implications to 
cosmological studies. Lilly et al. (1996) found the rapid rise of the comoving luminosity 
density with redshift in ultraviolet (UV), optical, and near-IR wavelenghts up to z ~ 1. 
Madau et al. (1996) interpreted this as an increase of comoving star formation rate (SFR) 
by a factor of ten in the redshift interval 2 = — 1. Subsequent observational evidences 
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support this conclusion (e.g. Ellis et al. 1996; Tresse & Maddox 1998). Because FIR is 
the wavelength at which dust reradiates the optical - UV radiation from young stars, the 
steeply increasing SFR density is expected to propose a significant increment of FIR source 
counts, which may be related to the evolving population of star-forming galaxies discovered 
by Infrared Astronomical Satellite (IRAS) (Hacking & Houck 1987 (HH87); Hacking, 
Condon, & Houck 1987; Saunders et al. 1990; Oliver, Rowan-Robinson, & Saunders 1992; 
Ashby et al. 1996). 

For the understanding of FIR properties of actively star-forming galaxies, vast advances 
have been provided by the IRAS all-sky survey with a flux limit of ~ 1 Jy at 60 /im (for a 
review, Soifer, Houck, & Neugebauer 1987a). Recently quite important results have been 
reported from the Infrared Space Observatory (ISO; Kessler et al. 1996). Kawara et al. 
1998 (K98) have performed a deep survey at the Lockman Hole region with a flux limit of 
45 mJy at 175 /im. They reported that the surface density of the sources brighter than 
150 mJy at 175 /im approximately reproduced by the model of Guiderdoni et al. (1997). 
Although this survey is much deeper than the IRAS survey, it covers only a small solid 
angle of the sky (1600 arcmin 2 ) and much wider sky coverage is necessary in the next step 
to obtain a huge sample. 

Infrared Imaging Surveyor (ffiJSQ) is a next generation infrared satellite, which will be 
launched in the beginning of 2003. A near-/mid-infrared camera (IRC) and a FIR scanner 
(FIS) are planned to be on-board. One of the main purposes of the IRIS mission is a FIR 
all-sky survey deeper than IRAS surveys, and the FIR instrument represents a significant 
improvement over IRAS, so that much higher sensitivity is achieved (detection limit is ~ 20 
mJy at 50 /mi, ~ 30 mJy at 120 /mi, and ~ 50 mJy at 150 /im). 

In order to examine the performance of the survey, we estimated the FIR galaxy counts 
based on a simple model with various sets of cosmological parameters and evolution types. 
There have been a number of attempts to predict the source counts or the CIRB properties. 
The models of the source counts and CIRB generally fall into two categories (e.g. Eales 
& Edmunds 1997): one is tied closely to the nearby observational results but makes no 
attempts to incorporate detailed galaxy evolution models (e.g. Beichman & Helou 1991 
(BH91); Pearson & Rowan- Robinson 1996). These models are based on the local FIR 
luminosity function (LF) and the observational evidence of spectral energy distribution 
(SED) of galaxies in the FIR - sub-mm, with the assumptions of simple functional forms 
for the evolution. The other is based on the models constructed by detailed physical 
processes related to the evolution of galaxies (e.g. Franceschini et al. 1994; Guiderdoni et 
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al. 1997). They tend to be very complex with large numbers of parameters, which give 
rise to a number of related problems poorly understood. Ellis (1997) names the former 
"empirical approach" and the latter "ab initio approach" . Despite each category of models 
has its own merits and demerits, we preferred the "empirical approach", because it is free 
of the parameters which are inherent in specific models and thus easy to understand the 
dependence on the cosmological and evolutional parameters. 

This paper is organized as follows. In section 2, we describe the general formulation of 
our model. The parameters to execute our calculations and their consistency with previous 
survey results are examined in section 3. We present the JRZS-based results and related 
discussions in section 4, and our conclusions are summarized in section 5. 



2. MODEL DESCRIPTION 

2.1. Spectral Energy Distribution of Galaxies 

First we construct a model of galaxy SED from FIR to sub-mm wavelengths. The FIR 
emission originates from several dust components of different grain sizes and temperatures, 
reprocessing the interstellar radiation (Desert, Boulanger, & Puget 1990). Rowan- Robinson 
& Crawford (1989) (RC89) showed that the FIR spectrum of IRAS galaxies is well fitted 
by the combination of two components: cool disk emission from interstellar dust (cirrus 
component) and hot starburst-induced emission, from optically thick dust clouds heated by 
numbers of OB stars. Therefore we consider a simple model of two components, cool cirrus 
and hot starbursts as the components of our model SEDs of galaxies. We must note that 
we ignore the power-law component seen in AGN or QSO spectra, because contribution of 
these power-law objects is expected to be negligible in number compared with the total 
number of infrared (IR) sources, and because their luminosity contribution is also small in 
FIR wavelength range, though it might be slightly higher in very high-z Universe (Pearson 
1996). 

We define the total IR luminosity L m as follows: 

/"max 
L(v)dv, (1) 
- 'min 

where L(u) is the monochromatic luminosity of a galaxy, and we put iVm = 3 x 10 11 Hz 
(A max = 1 mm) and z/ max = 10 14 Hz (A min = 3 //m). 

Here we assume that the IR luminosity of a galaxy less than 1O 1O L is composed of the 
cirrus component L c (u) only. For galaxies with L IR > 10 10 L Q , the IR luminosity in excess 
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of 10 10 L Q (i.e., L IR — 1O 1O L ), is assumed to come from the starburst component, denoted 
as L s (v). Thus, the composite SED of a galaxy is expressed as follows; 

T( 1 A = \ L ^ for L IR < 1O 1O L ; 

1 ) \ L c (u) + L s (u) for L m > 1O 1O L , 1 ) 



where LJv) and L s (u) are normalized as 



, _ / r / i i J ^ IR for Lir < 1O 1O L ; 

Lc= / ^("^H i O io L0 for Lir > 1O 1O L , (3) 



and 



- T ( , A JO for L m < 1O 1U L ; 

n Ls(Z/)dZ/ = lL IR -lO-L for L IR > 1O-L . (4) 

This model formulation is in line of BH91. So as to make our model work well, Lpr should 
be uniquely defined from the value L(v) at any frequency v. The SEDs depend only on 
Lju, and do not vary with redshift as long as of a galaxy is constant. 



2.1.1. The cirrus component 

The cirrus spectrum we used is based on the model of Desert et al. (1990). Their model 
spectrum is the sum of three subcomponents, that are polycyclic aromatic hydrocarbons 
(PAHs; peaked at ~ 10 /mi), very small grains (peaked at ~ 60 /im), and big grains (peaked 
at ~ 100 /im). No IR emission lines are included in our model. As mentioned above, 
galaxies with < 1O 1O L only have the cirrus component in the present model. At each 
IR luminosity, the galaxy SED is scaled properly to yield the given Lir according to eq.(|3]). 



2.1.2. The Starburst Component 

As for galaxies with L IR ^ 1O 1O ~ U L , bulk of the IR luminosity is due to dust 
emission from an intense starburst in giant molecular clouds (e.g. Sanders & Mirabel 
1996). Therefore we assume that, for galaxies whose IR luminosities are brighter than 
1O 1O L , their IR luminosities in excess of 1O 1O L consist of starbursts with two different 
temperatures according to BH91. We give the starburst flux L s (u) as a superposition of 
thermal blackbody radiation spectrum of temperature Th ot (we call it "hot component" 
of starbursts) and that of T coo i ("cool component"), multiplied with dust emissivity. We 
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adopted the wavelength-dependent dust emissivity oc A 1 oc v (RC89). By using the above 
temperatures, L a {y) is expressed by 

L s (u) = avB v (T cool ) + f3vB v (T hot ), (5) 

where B U (T) is the Planck function with temperature T and both a and f3 are normalizing 
constants (RC89). The temperatures of the starburst components are given by 

( L \ 01 

Tcooi = 60 — -V- K, (6) 



lO 11 ^ 



0, 



and 



/ L 



o.i 



(BH91; see also Rieke & Lebofsky 1986; Helou 1986; Soifer et al. 1987b). The cool 
component gives a good representation of the larger (0.01-0.1 /im) grains modeled by 
Rowan- Robinson (1986), and the hot component gives an approximate representation of the 
very small grains postulated by Boulanger, Baud, & van Albada (1985), Draine & Anderson 
(1985), and Rowan-Robinson (1992). 

Normalizing constants, a and (3, are determined by the following equations (BH91): 

/•^max 

0.7L S = a uB u (T cool ) du, (8) 

J ^min 

0.3L S = / isB u (T hot )dis. (9) 



The model SEDs of galaxies with Lir = 1O 8 L — 1O 14 L constructed in this manner 
are shown in Fig. 1. The solid curves represent our model SEDs with various Lir 
(the lowest one: Lfr = 1O 8 L ; level interval: AlogLi R = 1.0). We also presented the 
observed SEDs of three representative IR galaxies, M82 (Klein, Wielebinski, & Morsi 1988), 
IRAS F10214 + 4724 (z = 2.286, Rowan-Robinson et al. 1993), and SMM02399 - 0136 
(z = 2.803, Ivison et al. 1998). Open triangles stand for the observed SED of M82, open 
squares depict the SED of IRAS F10214 + 4724, and open diamonds represent the SED of 
SMM02399 — 0136. Small downward arrows show upper-limit values. The IR luminosity 
of M82 is ~ 4 x 1O 1O L and that of IRAS F10214 + 4724 is ~ 3 x 1O 14 L , but the latter 
has been revealed to be magnified by gravitational lensing by a factor of ~ 20 (Broadhurst 
& Lehar 1995). Its L m turns out to be ~ 1O 13 L , and its SED is corrected in Fig. 1. 
SMM02399 — 0136 is also a lensed object, and its unlensed L m is ~ 4 x 10 13 (Ivison et al. 
1998). The SED of this galaxy is also corrected for the lensing (by factors of 3.5). Their 
SEDs are approximately fitted by our models with corresponding L IR . 
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2.2. Local Luminosity Function 

We adopted the LF based on the IRAS by Soifer et al. (1987b) as the local IR LF of 
galaxies . Their LF is defined by L e0fim , so we converted it into our L IR , using our model 
SEDs. Then we performed an analytical fitting to the resultant data points, and obtained 
the following double power-law form for our LF: 

f 7.9 - 1.0 log (h^L m /L e ) for 1O 8 L < h^L m < 1O 1O - 3 L ; 
\og[h 3 75 Mhr 5 2 Lm)} = 17.1 - 1.9 log (h^L lK /L Q ) for 1O 1O - 3 L < h^L m < 1O 14 L ; (10) 

[ no galaxies otherwise, 

where 0o is the number density of galaxies in Mpc~ 3 dex" 1 in case of = 1, where h 75 is 
the Hubble parameter normalized by 75 km s~ x Mpc -1 (i.e. ^75 = Hq/75 km s _1 Mpc -1 ). 
This LF is presented in Fig. 2. We note that, because the volume element (defined in eq.^f]) 
depends on H^ 3 , luminosity function multiplied by the volume element has no dependence 
on the Hubble parameter. This means that neither galaxy number count nor the expected 
redshift distribution depend on the Hubble parameter in the empirical approach. Thus, we 
do not have to consider the effect of the Hubble parameter in the following analyses. 

The faint-end slope of the LF of Soifer et al. (1987b) is much steeper than that derived 
by Saunders et al. (1990). Little is known about the population which occupies the faint 
end of the LF. We discuss the effect of it in sections |2.5| and [4.3| . 



2.3. Treatment of Galaxy Evolution 

For galaxy evolution, we deal with the following three cases: (i) no evolution (ii) pure 
luminosity evolution, and (iii) pure density evolution. We define the LF at z, (p(z, Lir) in 
comoving volume. Then various types of galaxy evolution are characterized by the form of 
the functions which represent the relation between O (L IR ) and <p(z , Lin). We also define 
cumulative LF at redshift z, $(2, -^ir), 



L m ) = / 4>(z, L' lR ) d log Lj R (11) 
= / 4>(z, £rn)dlogL IR , 

J Lib. 



where LiR max is the maximum IR luminosity at the epoch. 
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2.3.1. No evolution 

The term "no evolution" means that the luminosity function is independent of z: 

4>{z, L m ) = MLm)- (12) 
Thus the cumulative LF reduces to 

PCX) 

*(z,Lir)=/ ML' m )d\ogL' m . (13) 

Although no evolution model seems to describe an unphysical situation, its prediction 
acts as a valuable standard baseline from which the various evolutionary differences can be 
compared. 



2.3.2. Pure luminosity evolution 



The effect of the luminosity evolution of galaxies is modeled by 

L m (z) = L m (0)f(z) , 



f(z) = exp 



Q 



r(*) 
tn 



(14) 
(15) 



where Q, t(z) and tn are a parameter defining the magnitude of evolution, look-back time 
as a function of z, and the Hubble time l/H . This functional form of evolution is proposed 
by Broadhurst, Ellis, & Glazebrook (1992). Peacock (1987) derived a useful approximation 
for the look-back time: 

T(Z) + (16) 



t H P 



where 



P = l + 



(2g ) 



0.6 



(17) 



This expression for the look-back time is exact for q = and q = 1/2 and is accurate 
within the error of < 1 % for q < 3/2 (Peacock 1987). The function reduces to the simpler 
form f(z) ~ (1 + z) Q at the \ow-z epoch. 

Pure luminosity evolution leads to the following expressions as cumulative LF, 

roo 

$(z,L m ) = / (f)(z,L' m ) d\ogL' m 

J Ltr 

^(L' m )d\ogL' m , (18) 



^IR 



Lmjf(z) 



(19) 
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2.3.3. Pure density evolution 

We, here, consider the model that only the comoving number density of galaxies 
chang GS clS £L function of redshift, which is called "pure density evolution". Pure density 
evolution is expressed as 

<P(z,L m )=MLm)g(z) ■ (20) 

The function g(z) is defined by 

(21) 

where P is a parameter defining the amplitude of evolution. This is the same type of 
function considered in the previous subsection. We note that this type of galaxy evolution 
is not aimed at describing the "merging" process. The number count model including 
merging process is widely discussed by Gardner (1998). The cumulative LF becomes 

<S>(z, L m ) = / ML' 1R )g(z)d\ogL' m . (22) 

J Lin 



g(z) = exp 



2.4. Galaxy Number Count 

Using the above formulae, now we calculate the flux-number relation, or so-called 
galaxy number count. We assume that galaxies are regarded as point sources (i.e. 
cosmological dimming of surface brightness is not taken into account). Then the relation 
between observed flux S(v) and emitted monochromatic luminosity L{u) is given by 



[l + z)L((l_+z)u) 

11 



where is luminosity distance. In the Universe without the cosmological constant A, g?l is 
expressed by analytical form, in terms of redshift z, deceleration parameter go and Hubble 
parameter H (Mattig 1958); 



d L = -^[zq + (g - l)(v/2g z + l- 1)] , (24) 

c : the velocity of light). When we fix a certain S(v), we obtain L((l + z)v) by using eq. 
|H). Then the correspondent L IR (S'(z/), z) at the redshift z is uniquely determined (see 
TU 
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We define N(> S(u)) as the number of galaxies with a detected flux density larger 
than S(u), then it is formulated as 

r f2 max d^V t°° 

N(> S(u)) = I ' dfl / dz— - 6(z, L' m ) d\ogL' m , (25) 

Jn Jo az ds I Jl 1k (S(l>) ,z) 

where d 2 ^/dzdf2 is the comoving volume element per str per z, which can be expressed in 
terms of cosmological parameters as 



d 2 V 



dttdz 



H (1 + zfy/l + 2q z 



dl (26) 



(e.g. Kolb & Turner 1994). We set z max as 10. 

Next, we formulate the number-flux relation with our evolution models. The 
no-evolution prediction is as follows: 

N(> 5(1/)) = / dO / dz— - / ML' m ) dlogL; R . (27) 

Jn Jo az a\i Jl ir (s(u) ,z) 

If pure luminosity evolution takes place, the expected number count is derived from eqs.(|i~8|) 
and (|D , 

N(> SM) = / dtt / d^— — / ^o(^r) dlogL' . (28) 

v K JJ Jn Jo dzdnJ Lm (S(u),z)/m y v m; & ir y 1 

The expression for pure density evolution is obtained in the same manner; 

f /"Zmax d?V I" 00 

N(> S{u)) = / dO / dz—- <f> (L' lR )g(z) d\o S L' m . (29) 

Jn Jo dz diz JL m (s(v),z) 



2.5. Redshift Distribution 

When the flux detection limit Sn m (i>) is given, the limiting luminosity at redshift z 
can be obtained by the following process. From equation (|23|), the limiting monochromatic 
luminosity Li im (u em ), where z/ em is the emitted frequency at the rest frame of a galaxy at 
redshift z, can be calculated as 

-Miml^emJ — ^ ~ Jlim I T/T~^ J ' ^ / 

We, then, obtain the total IR luminosity Lm,ii m , which uniquely corresponds to the given 
-kiim^em), and galaxies with Ljr larger than Lir u m are detected. In terms of the derived 
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Lir lim, the redshift distribution is formulated as 



d 2 N 



dttdz 



d 2 V 



dn dz 
d 2 V 



{z,L' m )d\ogL' lR 



dn dz 



$(z,L 



IR, lim J 



(31) 



The evolutionary effect can be easily introduced in the redshift distribution by using eqs. 
( [14]) or (p0|), as <&(z, Lir) in equation (|3~TD, and substituting LiR 5 ii m instead of Lir. 



3. PARAMETERS FOR THE EVOLUTIONARY MODELS 

In this section, we estimate the evolutionary parameters P and Q defined in section 
|2.3| in order to perform our calculation. Then we check the validity and consistency of our 
predictions by comparison with the IRAS- and ISO-based number counts and the CIRB 
spectrum obtained by COBE . 



3.1. Evolutionary Parameters 

We used IRAS extragalactic source count data to obtain the values of parameters P 
and Q. Among IRAS bandpasses (12, 25, 60, and 100 /im), the 60- /im band is known to 
be the most suitable for extragalactic studies, and intensively used in the context of the 
investigation of galaxies (e.g. Soifer et al. 1987a). Therefore we employed the 60- /zm data 
of IRAS galaxies. 

We first present the differential flux-number relation of IRAS 60-/iin databases, 
IRAS Point Source Catalog (IRAS PSC) and HH87. This is depicted in Fig. 3. Open 
squares represent the data points of IRAS PSC, and filled squares, HH87. These diagrams 
are constant against flux if the Universe is static Euclidean and source distribution 
is homogeneous. The prediction of no-evolution model is also shown in Fig. 3. We 
perform least-square fitting of our evolutionary models on these data points. When 
we use all data points with equal weight, the most likely values for the parameters are 
(P, Q) = (1.3, 0.7) (g = 0.1) and (1.4, 0.8) (g = 0.5). We obtain higher values 
(P, Q) = (2.7, 1.4) (g = 0.1) and (2.8, 1.5) (g = 0.5) when the four statistically poorest 
points around logS [Jy] ~ —0.5 are omitted. We adopt the latter for our studies on 
evolution. We note that the estimation of P and Q depends on the fiducial point which is 
used to determine the normalization of galaxy counts. Some authors define the parameters 
P and Q so that the model predictions have the same value of (diV/ dS) S 2 ' 5 at 1 Jy (see 
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e.g., Fig. 1 of Oliver et al. 1992). The others estimated P and Q by overall fitting, as the 
same as we did, and did not made such kind of normalization (see e.g., Fig. 5 of Ashby 
et al. 1996). Therefore our values for P and Q seem smaller than those who adopted the 
former method, but the effect of the evolution is almost the same. These parameters are 
summarized in Table 1. They are almost independent of qo, at least within the available 
flux limit of IRAS data. 



3.2. Galaxy Count : Comparison with IRAS and ISO Results 

We compared our number count predictions with the previous survey results obtained 
by IRAS and ISO. In the upper two panels of Fig. 4, the 60-yum number count predictions 
with two evolutionary models (go = 0.1 and 0.5) are presented with the result of QMW 
IRAS galaxy survey (Rowan- Robinson et al. 1991). Within the depth of the QMW catalog 
there is only a small difference between no-evolution and evolutionary predictions, and the 
model predictions agree with the observed counts. The evolutionary effect is significant in 
fainter flux regime (e.g. Bertin, Dennefeld, & Moshir 1997). 

Comparison between our 175-/im galaxy number count and the results of the 
ISO Lockman Hole survey (K98) is illustrated in the lower panels of Fig. 4. Source 
identification and flux calibration contains considerable difficulties, which turned out to be 
an uncertainty of the number count by a few factors. Both IRAS and ISO number counts 
are successfully reproduced by our model calculations. 



3.3. Cosmic Infrared Background 

The CIRB is generated from the integrated light of galaxies. Therefore, combining 
the SEDs of galaxies and number count predictions, we obtain the CIRB spectrum. The 
observed flux density of a galaxy whose IR luminosity is Lm, S(u, Lir) is given by equation 
( p3|) as follows: 

q(v r v (l + z)L(u(l+z),L m ) 

S{v , L m ) ^ , (32) 

where L(v, Lm) is the monochromatic luminosity of a galaxy with Lm. Then the CIRB 
spectrum 1(f), i.e., the background flux density from unit solid angle, is expressed as 

i(y) = I dz a^n L ^ L '™ )s{u > dlogL '™ • (33) 
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We can deal with the evolutionary effect on I(y) through (p(z, L IR ), just the same 
as in the case of number count calculation. The expected CIRB spectra are presented 
in Fig. 5. The upper panel shows the case of go — 0.1, and lower panel, = 0.5. We 
also put the observational constraints on I{y) obtained by COBE measurement in Fig. 
5. Red open triangles represent the "dark sky" upper limit to the CIRB measured by 
COBE Diffuse Infrared Background Experiment (DIRBE). The DIRBE sky brightness 
varies roughly sinusoidally over the year, due to the complex features of the interplanetary 
dust cloud. The "dark sky" brightness is that of the darkest area on the sky at each 
wavelength. Filled squares are the result of a similar analysis of the COBE Far Infrared 
Absolute Spectrophotometer (FIRAS) high-frequency data, after removal of the CMB 
(cosmic microwave background) signal. The dark sky values of DIRBE and FIRAS show 
excellent agreement with each other. Blue open triangles represent the residual signal 
of DIRBE after removing the contributions from the model foreground sources such as 
Galactic diffuse emission, and interplanetary dust emission. As described by Mather et 
al.(1994), the CMB spectrum in the wavelength range 0.5 — 5 mm deviates from a 2.726 
K-blackbody shape by less than 0.03 % of the peak intensity. This deviation is shown 
by green horizontal lines. Details of these background signals are summarized in Hauser 
(1995). In addition, Puget et al. (1996) have reported the isotropic diffuse residual light 
of FIRAS data, after subtracting the foreground contributions. Fixsen et al. (1998) and 
Hauser et al. (1998) have also reported similar results. The orange broken solid lines are 
the approximate shape of the spectrum of the residual light reported by Puget et al., Fixsen 
et al., and Hauser et al. with errors. 

Within the wavelengths shorter than 300 /im, our model prediction is consistent with 
the CIRB intensity, but in the longer wavelengths, it seems to violate the limits from 
the CIRB (larger q mitigates this discrepancy). In such a long wavelength regime, the 
background radiation is dominated by the contributions from hyperluminous IR galaxies 
at very high redshifts. As we see from Fig. 1, the SEDs of such luminous galaxies are 
well represented by our SED models. Thus, the functional form of the galaxy evolution 
might be too simple and possibly not valid in extremely high- z objects, or alternatively, 
high-luminosity end of the LF might not be valid at high-z. In spite of the difficulty in 
sub-mm range, our model works for our present purpose which is IR source number count 
prediction in wavelength region shorter than 300 /im. 

4. RESULTS AND DISCUSSIONS 
4.1. IRIS Survey Performance 
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4-1.1. Flux Detection Limit of IRIS FIS 

The flux detection limit Sn m (u) is decided by some factors; the sensitivity of the device, 
the intervening emission, e.g. zodiacal light, interstellar emission of Our Galaxy, and source 
confusion etc. The dominant factor is different in different wavebands, and the resultant 
S\im(v) is dependent on wavelength. Thus the calculation requires the careful consideration 
of S]\ m (u) individually for each A or v. We assume 50-, 70-, 120-, and 150/im bandpasses for 
IRIS FIS. Unstressed gallium-doped germanium (Ge:Ga) is used for the shorter wavelength 
in FIR (50 — 110 /mi) and stressed Ge:Ga for the longer wavelength (< 170 //m) (Kawada 
et al. 1998). In the shorter wavelength, sensitivity is limited by internal and background 
noises, while in the longer wavelength, it is constrained by the source and Galactic cirrus 
confusion (Gautier et al. 1992; Thronson et al. 1995). The resultant 5<r-detection limits of 
the assumed four wavebands are summarized in Table 2Q. 



4-1.2. IRIS galaxy number count predictions 

Our predictions for the galaxy number count at assumed IRIS four bandpasses are 
shown in Fig. 6. The upper four panels depict the number counts with go — 0.1 and 
the lower four, q = 0.5; from left to right, A = 50, 70, 120, and 150 /im. The black solid 
lines show the no-evolution prediction. The red and cyan lines represent the case of pure 
luminosity evolution and pure density evolution, respectively. The green dot-dashed lines 
in Fig. 6 denote the IRIS FIS flux detection limit at each band. The integrated detection 
number to the flux limit of each bandpass is ~ 10 5 str -1 at 50 /im, ~ 2 x 10 5 str -1 at 70 /im, 
~ 3 x 10 5 str -1 at 120 /mi, and ~ 2 x 10 5 str -1 at 150 /mi in the case of no evolution, thus 
the expected detection number is largest at 120 /im, and the total number of galaxy is 
expected to be ~ several x 10 6 in the whole sky. This is about 100-times larger than the 
number of IRAS galaxies, which is ~ 25000 in all the sky. It is clearly seen that go does not 
affect the number count within the detectable flux range of IRIS as shown in Fig. 6. 

The effect of evolution appears most significantly at 70 /mi among the assumed bands, 
and the total detection number becomes twice larger. Thus, the determination of the 
amplitude of galaxy evolution is best estimated at around 70 /im in the wavelength and 



2 The flux detection limits presented here are based on those in Kawada et al. (1998), but are slightly 
(by a factor of 1.5 - 2) deeper than the limits they reported; the values adopted here are those expected in 
the high Galactic latitude region (|6| > 60°). Much better point source detection limits can be expected in 
limited sky areas near the ecliptic poles, where the survey scan will be repeated more than hundred times, 
and can be achieved by a spatial deconvolution algorithm. 
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flux coverage of IRIS. With the aid of such a vast number of galaxy detections, we will 
be able to evaluate the evolution strength (P or Q) precisely. In the case of 70-/mi band, 
statistical fluctuation becomes much smaller than those of previous area-limited deep 
surveys. Consequently the errors in the estimation of P or Q will be small enough, at 
least, to determine them in the relatively low-redshift (z < 1) universe. Unfortunately, 
there exists no drastic difference between luminosity and density evolution in the behaviors 
of number count, and judging which evolution takes place from number count alone still 
remains as a difficult work. 

4.2. Redshift Distribution 

Using the IRIS FIS flux detection limits we gave in section |4.1.1] , we calculate the 
redshift distribution of galaxies detected by the IRIS all-sky survey. Figure 7 shows the 
expected redshift distribution of galaxies detected in the assumed four bands. Figure 7a 
illustrates the prediction at 50-/im band. Figures 7b, 7c, and 7d are the same as Fig. 
7a, except that they show the predictions at 70-, 120-, and 150-/mi bands, respectively. 
Solid curves denote the predictions for q = 0.1, and dashed curves, go = 0.5. Larger q 
leads to deeper ^-distribution. In any bandpasses, the numbers of galaxies at z > 1 are 
~ a few x 10 3 , and several times larger with evolution. Namely, low- ,2 galaxies are the 
overwhelming majority among the detected objects. The steep increases of galaxy number 
at z < 0.5 seen in Figs. 7c and 7d are attributed to the peak at ~ 100 /mi, which originates 
from big grains, in the SED of low-L IR (L IR <, 10 11 L ) galaxies. Figure 8 is an example 
of the redshift distribution for detected galaxies with various L IR (g = 0.1, no evolution). 
Such galaxies can be detected only within z ^ 0.5. The LF we employed makes them to be 
abundant; if we use flatter faint-end slope of the LF, the increases would be less prominent. 

The overall shape of the redshift distribution depends on the observing bandpass 
wavelength. Nearest galaxies (z ~ — 0.5) are more likely to be detected in longer 
wavelength A = 120 — 150 /mi, and their detection number steeply decreases with increasing 
z in the longest wavelengths. Galaxies which lie within the redshift range of z ~ 0.5 — 2.5 
are expected to be detected more at 70 /mi than at the other bands. When the galaxy 
redshift exceeds ~ 3, they become harder to be detected in the shorter bands, while 
they remain detectable in the longer bands. The detection numbers at A = 120 /mi and 
150 /mi, therefore, surpass those at A = 50 /mi and 70 /mi for the furthest galaxies. This 
is understood as follows. At significant redshift (z > 3), the detected objects are the most 
luminous starbursts, whose SEDs peak at around 30 /tm in the rest frame. Since the peaks 
of luminous starbursts shift away from their rest wavelength, they are not detected at 
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high-2; in shorter wavelengths. Instead, in the longer wavebands, redshift effect makes their 
observed flux density nearly constant (known as the negative iCcorrection), and they can 
be still seen at high-z. 



4.3. Color— color Diagram in FIR 



As we saw in section |4.2| , most of the objects detected in the IRIS galaxy survey are 
nearby ones. Hence, when we focus on the very young or primeval galaxies, an efficient 
method to pick up such candidates without knowing their spectroscopic redshifts will be 
very useful. We, here, make attempts to utilize the peaks of the SEDs to estimate the 
redshifts roughly from their FIR colors alone. We suppose that all the Galactic objects are 
taken away by certain methods. This may well be possible at high Galactic latitudes. 

Figure 9 shows the loci of galaxies with various L IR on the FIR color-color (logSgo/SVo 
- log S50/S150) plane, where S\ is the detected flux density of galaxies at wavelength A 
[^m] in unit of [erg s _1 cm _2 Hz~ 1 ]. Filled symbols are put on every Az = 0.5 interval from 
z = 0. When the sources fade away from the IRIS detection limits, we do not show their 
loci and symbols. Red filled squares represent the colors of the \ow-z (z < 1) galaxies, 
green filled squares the intermediate- z (1 < z < 2.5), and blue filled circles the high- z 
(2.5 < z) ones. Because of the way of the construction of our model SEDs, galaxies with 
Lir = 10 8 — 1O 1O L have the same colors (the red filled square at the top-right). To see 
how the existence of starburst component changes the color of cool cirrus galaxies, we put 
the galaxies of L IR = 1O 1O1 L — 10 10,9 L Q with an interval of AlogL IR = 0.1 (denoted by 
open squares from the top-right to down-left). They correspond to the galaxies whose IR 
luminosity fraction of the superposed starburst is 20.6, 36.7, 50.0, 60.2, 68.4, 74.9, 80.0, 
84.2, and 87.4 %, respectively. 

When we consider the color selection criteria, the intrinsic scatter of galaxy color has to 
be taken into account. Schmitt et al. (1997) proposed the template SEDs of various classes 
of galaxies. Their SEDs show that the dispersion at FIR wavelength is an order unity. The 
scatter of IRAS Sioo/Seo is also known to be an order of magnitude (RC89). Considering 
these scatters, we can estimate the domains where galaxies locate on the FIR color-color 
plane. We show the the domains of low- (z < 1), intermediate- (1 < z < 2.5) and high-z 
(z > 2.5) galaxies with a red dashed rhomboid, a green solid rectangle, and a blue solid 
rectangle, respectively. The high-z and mid- z regions are well separated on this diagram. 

We must think of the overlap between low- z and mid/high- z domains. The 
low- z population creeping in the mid/high- ,2 domain consists of only low-L IR galaxies 
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(L IR <, 10 11 L Q ). As we have seen in section such galaxies quickly become undetectable 
with increasing redshift and reside in z < 0.2. The expected optical magnitude based on 
SEDs extended to optical region ranges from ~ 17 mag to ~ 19 mag (Hirashita et al. 1998). 
Since the positional accuracy of the IRIS survey is quite good (5" - 10") and surface density 
of such bright galaxies is low enough, these galaxies would be identified easily with optical 
counterparts. On the contrary, the high-z galaxies are expected to have fainter magnitudes 
(~ 20 mag) (Hirashita et al. 1998). Thus the cross correlation with optical survey data such 
as POSS, or forthcoming SDSS would identify most of the low- 2: IR-galaxies. Evolutionary 
effect makes the mid/high- z galaxy more abundant, and the fraction of such objects in 
the overlapped region will be higher. Thus, color classification of galaxy redshift can be 
practicable using sets of bandpasses 50, 70, and 150 /im with the aid of present and future 
optical surveys. 

5. SUMMARY AND CONCLUSIONS 

In this paper we examined the performance of the IRIS FIR all-sky survey, one of whose 
main purposes is to study galaxy evolution and formation. We used a simple empirical 
model for galaxy number count estimation, by adopting a multicomponent SED model 
which consists of cirrus and starburst components, and the nearby FIR luminosity function 
derived from that of IRAS galaxies. 

Our conclusions are as follows: 

1. A few xlO 6 galaxies will be detected by the IRIS all-sky survey, which is 100 times 
larger number of detection than that of IRAS. 

2. Shorter wavelength bandpass (50 — 70 /im) is suitable for intermediate- z (1 < z < 2.5) 
galaxy detection, while longer wavelength (120 — 150 jum) one is for high-z (2.5 < z) 
and nearby (z < 1) galaxies. 

3. We will detect the effect of galaxy evolution and can evaluate the amplitude of 
evolution at least in the nearby universe in the IRIS survey, but still it is difficult 
to constrain which type of evolution takes place, from the number count alone. On 
the contrary, the intensity of the CIRB complementarily provides us much tighter 
constraint on the form of galaxy evolution at very high- z. 

4. Redshift estimation using galaxy FIR colors can be practicable using sets of 
bandpasses 50, 70, and 150 /im with the aid of present and future optical surveys. 
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A. APPENDIX : OTHER FUTURE MISSIONS 

We show the expected number counts and redshift distributions of galaxies in other 
future missions now planned, SOFIA, SIRTF, FIRST, and LMSA. The present model is 
valid only within the mid-IR to sub-mm wavelengths because of the SED used, we restrict 
the count predictions within this range. We calculated the number counts and redshift 
distributions in the same way as we did in the main text. In this appendix, we only consider 
the case of go — 0.1 cosmology with pure luminosity evolution (Q = 1.4). Though some 
of the facilities are not planned to be used for large-area survey, we uniformly calculated 
the number counts per unit solid angle, and the redshift distributions in the whole sky. We 
applied the detection limits from published papers or from the web pages of the missions if 
available. 



A.l. SOFIA 

Stratospheric Observatory for Infrared Astronomy (SOFIA) is an airborne observatory 
with a 2.5-meter telescope installed in a Boeing 747 aircraft, which will be flying in 2001. 
The detection limits (la noise level) for 1 hour exposure are 7 mJy at 100 /im and 3.5 
mJy at 450 //m (Becklin 1997). Based on these values we calculated the number count and 
redshifts with 5<r-detection limits (obtained by simply multiplying factor 5 to the above 
limits). The galaxy number counts are presented in the upper-left panel of Fig. 10, and the 
redshift distributions are shown in the upper-left panel of Fig. 11. 
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A.2. SIRTF 

The Space Infrared Telescope Facility (SIRTF) is the 4th and final element in NASA's 
family of "Great Observatories". SIRTF consists of a 0.85-meter telescope and three cooled 
instruments, IRAC, MIPS, and IRS, capable of performing imaging and spectroscopy in the 
3 — 180 (Km wavelength range. It is planned to be launched in December 2001. We present 
the predictions for MIPS (Multiband Imaging Photometer for SIRTF) at wavelengths 23.5, 
70, and 160 jum. The 5cr-detection limits are 370 /iJy, 1.4 mJy, and 7.5 mJy, respectively, 
by 500-s exposure (Heim et al. 1998). SIRTF will survey 60 % of the whole sky. The 
upper-left panel of Fig. 10 depicts the number counts, and the upper-left panel of Fig. 11 
shows the redshift distribution. 



A.3. FIRST 

The Far Infrared and Submillimetre Telescope (FIRST) is the 4th cornerstone mission 
in the European Space Agency's Horizons 2000 programme, implemented in collaboration 
with NASA. FIRST will perform photometry and spectroscopy in the 80 — 670 /zm range. 
Three instruments have been provisionally selected: HIFI, PACS, and SPIRE. Imaging 
photometry is performed by PACS (80 — 210 fim) and SPIRE (200 — 670 pm). The 
5cr-detection limits (by 1 hour exposure) for PACS and SPIRE are 5 mJy and 3 mJy, 
respectively (Pilbratt 1998). We show 90- /xm number count for PACS and 250, 350, and 
500 jum for SPIRE in the lower-left panel of Fig. 10. Redshift distributions of the objects 
detedted at these bandpasses are presented in the lower-left panel of Fig. 11. 

A.4. LMSA 

The Large Millimeter and Sub millimeter Array (LMSA) is the ground-based radio 
facility proposed in Japan. In the current design concept the array will consist of 50 
10-m antennas and will be covering observing frequencies from 80 to 800 GHz. The array 
will be located at very high site to realize sub-arcsec resolution imaging at very high 
frequencies. We present the number counts expected by LMSA at 350, 450, 650, and 850 
/im. The 5<r sensitivities at these wavelengths (yearly mean values) are 3350, 1800, 750, 
and 115 /iJy/beam, respectively, by 8- hour integration (Kawabe & Kohno 1998, private 
communication). Under the best weather condition in winter, much higher sensitivities can 
be achieved. The galaxy number counts are shown in the lower- right panel of Fig. 10. The 
lower-right panel of Fig. 11 illustrates the expected redshifts. In the lower-right panel of 
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Fig. 10, we compared our prediction with the 850 //m number counts recently obtained by 
SCUBA in JCMT. Filled red circle represents the value obtained by Smail et al. (1998), 
filled square shows the result in Hubble Deep Field (Hughes et al. 1998), and filled triangle 
is the number count in Lockman Hole (Barger et al. 1998). In spite that we include the 
evolution in this calculation, our prediction is an order of magnitude smaller than those 
results. Thus, a very strong evolution is suggested by sub-mm observations. It requires 
further detailed analyses of the present observational data to study the history of galaxy 
evolution in wide range of wavelength. 
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Figure Captions 



Fig. 1.— The model spectral energy distributions (SEDs) of galaxies with Lir = 
10 8 L Q — 1O 14 L . The solid curves represent our model SEDs with various Lir (the lowest 
one: Lir = 1O 8 L ; level interval: A log Lir = 1-0). We also present the observed SEDs of 
two representative IR galaxies, M82 (Klein, Wielebinski, & Morsi 1988), IRAS F10214+4724 
(Rowan-Robinson et al. 1993), and SMM02399 - 0136 (Ivison et al. 1998). Open triangles 
stand for the observed SED of M82, open squares depict the SED of IRAS F10214 + 4724, 
and open diamonds represent the SED of SMM02399 — 0136. Small downward arrows show 
upper-limit values. The IR luminosity of M82 is ~ 4 x 1O 1O L and that of IRAS F10214+4724 
is ~ 3 x 1O 14 L , but the latter has been revealed to be magnified by gravitational lensing by 
a factor of ~ 20 (Broadhurst & Lehar 1995). Its Lir turns out to be ~ 1O 13 L , and its SED 
is corrected in this figure. SMM02399 — 0136 is also a lensed object, and its unlensed Lir 
is ~ 4 x 10 13 (Ivison et al. 1998). The SED of this galaxy is also corrected for the lensing. 
Their SEDs are approximately fitted by our models with corresponding Lir. 

Fig. 2. — The local infrared luminosity function (LF) of galaxies 0o(-^ir)- This is constructed 
from the LF of Soifer et al. (1987b). Since their LF is defined by L 60Atm , we converted it into 
Lir (the IR luminosity in the wavelength range of 3 /im — 1 mm), using our model SEDs. 
The filled squares represent the converted data points, and solid line stands for the fitted 



analytical function presented in the main text (see eq. 10). Dotted line is shown for the 
Lir of the 'knee' of our LF. 

Fig. 3. — The differential flux-number relation for IRAS 60-yum databases and our model 
predictions. Open squares represent the data points for IRAS PSC, and filled squares, 
HH87. The no-evolution model tracks are also shown by black solid curves. Upper 
panel shows the case of g = 0.1, and lower panel, of q = 0.5. We performed least- 
square fitting of our evolutionary models (see section |2~3| ) on these data points. When 
we adopt all data points with equal weight, the most likely value for the parameters are 
(P, Q) = (1.3, 0.7) (g = 0.1) and (1.4, 0.8) (q = 0.5). We obtain higher values 
(P, Q) = (2.7, 1.4) (g = 0.1) and (2.8, 1.5) (go = 0.5) when the four statistically poorest 
points around logS [Jy] ~ —0.5 are removed. We adopt the latter in our studies. 

Fig. 4. — The integrated galaxy number count (g = 0.1 and 0.5) at 60-/zm and 175-/im. The 
results of QMW IRAS galaxy survey (Rowan-Robinson et al. 1991) and the ISO Lockman 
Hole survey (Kawara et al. 1998) are also presented. The upper panels represents the 60-/xm 
prediction. The agreement of the model and the QMW count is excellent. The lower panels 
shows the 175-/xm number count. The horizontal errors are caused by the uncertainty of flux 
calibration. Within the error, the model agrees with the ISO result. 
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Fig. 5. — The expected cosmic infrared background (CIRB) spectra with g = 0.1 and 
q = 0.5. The background intensity depends on g , so that larger q leads to smaller I{v). 
We also put the observational constraints on CIRB obtained by COBE measurement. Open 
red triangles represent the "dark sky" (see section |3.3| ) upper limit to the CIRB measured by 
COBE Diffuse Infrared Background Experiment (DIRBE). Filled squares are the result of 
a similar analysis of the COBE Far Infrared Absolute Spectrophotometer (FIRAS) high- 
frequency data, after removal of the cosmic microwave background signal. Open blue 
triangles represent the residual signal of DIRBE after removing the contributions from the 
model foreground sources such as Galactic diffuse emission, and interplanetary dust emission. 
The green horizontal lines represent the upper limit derived from Mather et al. (1994), and 
solid orange broken lines represent the approximate value of the residual light reported by 
Puget et al. (1996), Fixsen et al. (1998) and Hauser et al. (1998) with errors. The CIRB 
constraint on the magnitude of galaxy evolution is tighter than that obtained from the 
number count. 

Fig. 6. — The galaxy number count predictions at assumed IRIS four bandpasses. The 
upper four panels depict the number counts with q = 0.1 and the lower four, q = 0.5; 
from left to right, A = 50, 70, 120, and 150 //m. The black solid lines show the no-evolution 
prediction. The red lines and cyan lines represent the case of pure luminosity evolution and 
pure density evolution, respectively. The green dot-dashed lines denote the IRIS FIS flux 
detection limit at each waveband. 

Fig. 7. — The expected redshift distribution of galaxies detected in the assumed four bands. 
Figure 7a illustrates the prediction at 50-/^m band. Figures 7b, 7c, and 7d are the same as 
Fig. 7a, except that they show the predictions at 70-, 120-, and 150-/zm bands, respectively. 
Solid curves denote the predictions for q = 0.1, and dashed curves, go = 0.5. Larger g 
leads to deeper ^-distribution. The steep increases of galaxy number at z < 0.5 seen in Figs. 
7c and 7d are attributed to the peak at ~ 100 /im, which originates from big grains, in the 
SED of low-L IR galaxies. 

Fig. 8. — An example of the redshift distribution of detected galaxies with various 
(go = 0.1, no evolution). Solid curves represent the galaxies with Iow-Lir, dashed curves the 
intermediate-LiR, and dot-dashed curves the high-LiR. 

Fig. 9.— The loci of galaxies with various Lir on the FIR color-color (log S50/SV0 — 
log S50/S150) plane. Filled squares are put on every Az = 0.5 interval. Red filled squares 
represent the colors of the low- z (z < 1) galaxies, green filled squares the intermediate- 
z (1 < z < 2.5), and blue filled squares the high-z (2.5 < z) ones. For the sources 
fade away from the IRIS detection limits, we do not show their loci and symbols. To 
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see how the starburst changes the color of cool cirrus galaxies, we put the galaxies of 
L IR = 10 1IX1 L Q — 1O 1O 9 L with an interval of AlogL IR = 0.1 (denoted by red open squares 
from the top- right to down-left). The domain which is occupied mainly by \ow-z, Iow-Ljr 
galaxies are shown roughly by red dotted rhomboid. Green rectangle shows the domain of 
intermediate- z ones, and blue rectangle shows the region where high-z galaxies to be located. 

Fig. 10.— Integrated number counts of galaxies expected for SOFIA, SIRTF, FIRST, 
LMSA. SOFIA (upper-left) : Blue and red solid lines represent our 120 /im and 450 /im 
predictions, respectively. The vertical blue and red dot-dashed lines show the 5<r-detection 
limit of SOFIA with 1-hour exposure. We adopted go — 0.1 with pure luminosity evolution 
(Q = 1.4) in this and the following figures. SIRTF (upper-right) : Integrated number counts 
of galaxies expected for SIRTF. Blue, green, and red solid lines represent our 25, 70, and 160 
/im predictions, respectively. The vertical blue, green, and red dot-dashed lines show the 5cr- 
detection limit of SIRTF MIPS with 500-s exposure. FIRST (lower-left) : Integrated number 
counts of galaxies expected for FIRST. Blue, green, yellow and red solid lines represent our 
90, 250, 350, and 500 /im predictions, respectively. The vertical blue, green, yellow and 
red dot-dashed lines show the 5<x-detection limit of FIRST PACS (90 /an) and SPIRE (250, 
350, and 500 /im) with 1-hour exposure. LMSA (lower-right) : Integrated number counts 
of galaxies expected for LMSA. Blue, green, yellow and red solid lines represent our 350, 
450, 650, and 800 /im predictions, respectively. The vertical blue, green, yellow and red dot- 
dashed lines show the 5<r-detection limit of LMSA with 8-hour exposure. Filled red circle 
represents the galaxy number count reported by Smail et al. (1998), filled square shows the 
result in Hubble Deep Field (Hughes et al. 1998), and filled triangle is the count in Lockman 
Hole (Barger et al. 1998), all obtained with SCUBA in JCMT at 850 /im . 

Fig. 11.— The redshift distribution of galaxies detected with SOFIA, SIRTF MIPS, FIRST, 
and LMSA assuming the whole-sky observation; upper-left : SOFIA, upper-right : SIRTF, 
lower-left : FIRST, and lower-right : LMSA. 
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Table 1 : Adopted Evolutionary Parameters 
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Table 2 : Adopted Flux Limits of IRIS FIS 



Wavelength 
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[mJy] 
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Figure 2 
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Figure 4 
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Figure 9 
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Figure 11 



